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ABSTRACT 

Residual star formation at late times in early-type galaxies and their progenitors must be suppressed 
in order to explain the population of red, passively evolving systems we see today. Likewise, residual 
or newly accreted reservoirs of molecular gas that are fuelling star formation must be destroyed. This 
suppression of star formation in early-type galaxies is now commonly attributed to AGN feedback 
wherein the reservoir of gas is heated and expelled during a phase of accretion onto the central 
supermassive black hole. However, direct observational evidence for a link between the destruction 
of this molecular gas and an AGN phase has been missing so far. We present new mm- wavelength 
observations from the IRAM 30m telescope of a sample of low redshift SDSS early-type galaxies 
currently undergoing this process of quenching of latc-time star formation. Our observations show that 
the disappearance of the molecular gas coincides within less than 100 Myr with the onset of accretion 
onto the black hole and is too rapid to be due to star formation alone. Since our sample galaxies 
are not associated to powerful quasar activity or radio jets, we conclude that low-luminosity AGN 
episodes are sufficient to suppress residual star formation in early-type galaxies. This 'suppression 
mode' of AGN feedback is very different from the 'truncation mode' linked to powerful quasar activity 
during early phases of galaxy formation. 

Subject headings: galaxies: elliptical and lenticular, cD, galaxies: evolution, galaxies: formation, 
galaxies: active 



1. INTRODUCTION 

The low redshift galaxy population can be divided 
into blue sta rforming and red passively evolving systems 
(jBaldrv et al.l 120041) . The red population is dominated 
by elliptical and lenticular (early- type) galaxies, and 
in order to explain the properties of these systems 
as we observe them, almost all residual star forma- 
tion in them must be suppressed. This suppression 
must take place via the heating or expulsion of the 
molecular gas which is the fuel for star formation 
(|Dekel fc Silkl 119861 : iBenson et aIll2003D and is thought 
to be driven by the energy input from accreting super- 
massive black holes at the galaxies' cente r s in a process 
called AGN feedback (ISilk fc Rees' '199 81: iBenson et all 
jDiJyIatteo et all 12005,: .Cattaneo et all 120051: 
[20051 iSchawinski erall [20061 : 
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iKavirai et al.l l2007a': 'Khalat van et al. 1'2007^. 

The stellar populations of early-type galaxies indicate 
that the most massive galaxies are the first to undergo 
this process of suppression, forming their stars earliest 
and on the short e st time scales ([Thomas et al. 20051 
Nelan et a L| 120051; [ Bcrnardi et al.' 2006; iJimenez et all 
2007: Kaviraj et all 2008a; Kavira} .2008bl ). These in- 



tense episodes of star formation are followed by passive 
evolution with almost no subsequent star formation, al- 
though evoluti on via processes such as dry r nergers is still 
possible (e.g. iKhochfar fc Burkertl 120031 ; Ivan DokkumI 
|2005 (). It has been suggested that the mass at which 
galaxies undergo the transition from star-forming to 
quiescent de creases over cosmic time (jBell et al.l 120051 : 
iBundv et al.l [2006) , and so less massive galaxies have 
more extended histories of star formation ([Thomas et al.l 
[2005h . 

As new molecular gas is continuously supplied by 
accretion, mergers and stellar mass loss, most cur- 
rent models of galaxy formation suppress residual 
star formation by invoking a phase of accretion 
onto the supermassive black hole at the cen t er of 
the galaxy jSH k & Rees 1998; Cattan eo efall 120051 : 
ISchawinski et al.. .2006. : .Croton et al.i i20(j6[ ) as the en- 
ergy liberated by core-collapse supernoy ae is insufficient 
([Dekel fc Silklll986l : [Bmson et all [200 3'). The resulting 
active galactic nucleus (AGN) phase includes jets, radi- 
ation and outflows liberating sufficient energy to destroy 
the molecular gas reservoir by heating and expelling it, 
thus suppressing star formation by depriving the host 
galaxy of its fuel. Observations indicate that the mass 
of gal axies and their supermassive black h oles are corre - 
lated ([Ferrarese fc Merrittll2000l : iGebhardt et al.ll2000f) . 



2 



Schawinski et al. 



suggesting that the two co-evolve. The destruction of 
the gas reservoir may then in turn terminate the growth 
of the black hole as it deprives itself of material for fur- 
ther accretion as these systems transition to a main- 
tenance mode suppressing cooling f lows due to stellar 
mass loss and col d accretion (Ciotti &: Ostrike"^ Il997l . 
l2007t iBower et al.i r2006'). Alternative - or perhaps com- 
plementary - gas heating mechanisms to AGN feed- 
back have been prop o sed bylBirnboim et all (|2007D and 
iKhochfar fc Ostrikeil ()2008[ ). Investigating the role of 
AGN feedback directly is thus an important observa- 
tional goal. 

Early-type galaxies are known to pass through episodes 
of accretion onto their central black hole, but evidence 
which links this phase of their evolution to the removal 
of the molecular gas reservoir has until now been lacking. 
In this Paper, we present observational evidence for this 
process occurring in low- to intermediate-mass early-type 
galaxies at low redshift and identify low-luminosity AGN 
as the culprit. 

Throughout this work, we assume cosmological param- 
eters (rim = 0.3, Ha — 0.7, Hq — 70), consistent with the 
Wilkinson Microwave Anisotropy Probe (WMAP) Third 
Year results and their co mbination of results with other 
data (jSpergel et al.ll2007[ ). 

2. OBSERVING AGN FEEDBACK IN ACTION 

We aim to directly observe the effect of AGN on a 
galaxy's molecular gas as it changes fr om a starforming 
system to a passive early-type galaxy. iSchawinski et aLl 
([2007b), hereafter S07, have assembled a sample of ap- 
proximately 16,000 early-ty pe galaxies from the Sloan 
Digital Sky Survey (SDSS; lYork et all l2000l) selected 
by visual inspection in the redshift range of 0.05 < 
z < 0.10. Visual inspe ction of SP S S images is a 
powe r ful method (F uku gita et al.l |2004 iFukugita et al.l 
[2007t iLintott et all fBol and ensures the inclusion of 
early-type galaxies tha t are still experiencing star for- 
mation or host AGN (lYi et al.l l2005t ISchawinski et all 
I2007al: iKavirai et all l2007bf ). We identify the domi- 
nant sources of ionising radiation as star formation (SF), 
AGN, or both in each galaxy usin g ratios of nebular 
emission line s (iBaldwin et all 119811 : iMiller et al.l 120031 : 
iKewlev et al.ll2006D. We fit the stellar popula tions using 
models (jThomas et al.l 120031 : lMarastonll2005f ) and quan- 
tify the age and the mass-fraction of the last significant 
episode of star formation by combining broadband pho- 
tometry, from the UV to the near-infrared, with stellar 
absorption indices measured from optical SDSS spectra. 

The star formation history parameters used here are 
taken from S07, who use a two-burst model for the star 
formation history. The 'old' burst is modelled as a vari- 
able single stellar population (SSP) and the young burst 
as an exponentially declining burst with an e-folding time 
of 100 Myr. Tests show that even large shifts in the e- 
folding time do not generally change the best-fit young 
age and mass-fraction; only significantly shorter (~ 10 
Myr) or longer 1 Gyr) e-folding time shift the results 
and these produce worse fits to the data. The error bars 
in the young burst age ty (e.g. in Figure H]) include both 
formal errors from the fit to the photometry and spectral 
indices and errors introduced by the choice of parame- 
terization by marginalizing over the other parameters, 
including dust and metallicity. The marginalization over 



the variable 'old' age accounts for any variation of the 
star formation prior to the most recent event, which is 
accounted for by the young burst. 

S07 find that a large fraction of low- and intermediate- 
mass early-type galaxies (50 < a < 150 kms^^; 10^" < 
A^steiiar < 10^^ Mq) show & clear time sequence (S07; 
see also Figure |T]), which begins with actively starform- 
ing early- type galaxies, which formed typically a few per- 
cent of their stellar mass in the most recent episode of 
star formation. This is followed approximately 200 Myr 
later by a phase in which both AGN and star formation 
(AGN-I-SF) exist together. After the AGN-hSF phase, 
the luminosity of the AGN increases, while all traces of 
star formation disappear. The Seyfert AGN (SY) phase 
begins approximately 500 Myr after the beginning of the 
last episode of star formation. The AGN then declines 
in luminosity before the galaxy eventually settles into 
passive evolution in terms of its star formation history, 
though further evolution via dry mergers may still oc- 
cur. This striking coincidence between the aging of stel- 
lar populations and the rapid evolution of the nebular 
emission from being excited by star formation to being 
powered by black hole accretion suggests a role for the 
AGN in suppressing star formation. A possible caveat 
on the nature of the SF phase is the possibility that an 
AGN is present even in the first starforming phase, but is 
sufficiently obscured to leave the global ISM unaffected 
and thus does not influence the optical emission lines. 
Whether such an obscured AGN already suppresses (or 
perhaps even drives) star formation is unclear and further 
observations are necessary to investigate this question. 

Similar observation al results have also been presented 
bv lSalim et all (|2007f ). who discuss a possible evolution- 
ary sequence from star formation to quiesc ence via an 
AGN phase. The key difference between ISalim et al.l 
()2007| ) and S07 is that S07 consider specifically the early- 
type galaxy population and that S07 demonstrate that 
the galaxies along the sequence from star formation to 
quiescence via the AGN phas e form a genuine e volu- 
tionary sequence. Furthermore. I Graves et aP (12007! ) find 
that at a fixed velocity dispersion, early-type galaxies 
hosting LINERs are systematically younger than their 
quiescent counterparts, sampling the penultimate evolu- 
tionary phase of S07. The results of S07 favor an AGN 
feedback interpretation, but cannot rule out alternative 
explanations. 

If the AGN are responsible for the destruction of the 
molecular gas, this material should disappear abruptly 
with the appearance of the AGN in the AGN-I-SF and 
Seyfert phase. If instead the star formation in these blue 
early-type galaxies suffices to exhaust or heat the gas 
reservoir, we should see a steady decline of the amount 
of molecular gas with time. 

In order to test these competing hypotheses, we mea- 
sure the molecular gas content of a sample of galaxies 
along this time sequence. Molecular hydrogen is exceed- 
ingly difhcult to observe directly, so CO - which is an 
excellent tracer of molecular gas mass jSdo mon et al.l 
119921 : iNeininger et al.|[T998l ) - is commonly used instead. 

3. OBSERVATIONS 

From S07, we select early-type galaxies with stellar 
masses from 10^" — 10^^ Mq and randomly choose 10 
SF, 10 AGN+SF and 4 Seyfert AGN with velocity dis- 




stellar mass [M®] 




Fig. 1. — The color-mass relationship for the galaxies in our sample. We plot the galaxy stellar mass derived by spectral energy 
distribution fitting using the models of Maraston (20Q3), versus the optical u — r color. In each panel, morphological late-types are gray, 
quiescent early-types are orange and the various active early-types (classified by optical emission line ratios) are colored in in each panel. 
From left to right: blue (SF), green (AGN+SF), red (Seyfert AGN). We mar k the galaxies obser ved with the IRAM 30m telescope with 
larger points, number them and show their SDSS gri composite color images l ILupton et al.|[2004l) below in three blocks corresponding to 
the three classifications. This Figure illustrates the progression of early-type galaxies from the blue cloud of starforming galaxies via an 
AGN phase to the red sequence of passive galaxies. 



persion a < 120 km . The velocity dispersion limit en- 
sures that we select only targets from S07 that are of the 
type that take part in our AGN feedback time sequence. 
Objects with higher velocity dispersions typically have 
older stellar populations, and must have successfully sup- 
pressed cooling over long (Gyr) timescales via mecha- 
nisms such as LINERs, which are common in such pop- 
ulations. This sample of early-type galaxies spans the 
transformation from star formation via an AGN phase 
to quiescence and so provides us with the ideal labo- 
ratory to move from inferring to directly probing the 
role of AGN. The difference between selecting galaxies 
from within categories in this way and randomly select- 
ing galaxies from the entire sample is at a level which is 
much less than 2 sigma (assuming a binomial distribu- 
tion) and therefore does not produce a significant change 
in our results. 

We observed the C0(1 0) transition line in this sam- 
ple of galaxies along the time sequence with the IRAM 
(Institutio de Radioastronomfa Milimetrica) 30m tele- 
scope at Pico Veleta in July and December 2007. The 
IRAM 30M telescope allowed us to simultaneously also 
observe the C0(2 1) transition. We illustrate our 
sample selection on the color-mass diagram in Figure [T] 
and also provide SDSS images of our targets. The spec- 
tra were reduced with standard CLASS software and first 
order polynomial baselines are removed. In Figure [21 we 
show the reduced spectra, including the fits to line pro- 
files when the C0(1 — > 0) transition line is detected. We 
note that the gas masses derived from our observations 
are total gas masses. The size of the IRAM beam at the 
frequency of C0(1 ^ 0) is ~ 17.0", while the effective 



radii of our galaxies is on the order of a few arcseconds. 
We compute the CO luminosity of each target and con- 
vert it to a molecular gas mass using a conversio n con- 
stant of a = 1.5 K km s'^ pc^ (jEvans et al.ll2005D . The 
derived gas masses and further galaxy properties are in 
Table [T] The properties of the detected lines are listed 
in Tabled 

4. RESULTS 

The typical molecular gas masses of those nearby 
early-type galaxies with detections ha ve been esti i natecl 
to range between 10^ - 10^ Mq (jSage et alj l2007t 
ICombes et al.l 120071 ). corresponding to very low mass 
fractions. These previous studies focused on 'typi- 
cal' early- type galaxies draw n from the red sequence 
population (see also e .g. iFaber fc Gallagherl 119761 : 
iLake fc Schommeilll984l ) . While the gas cont ent of spiral 
AGN hosts has been studied in detail (e.g. iBlitz et al] 
Il986t iGarcia-Burillo et "all I2003D . no such studies exist 
for early-type galaxies. Our work is the first to system- 
atically observe gas in early-types systems by mapping 
the transition from blue to red via an AGN phase and 
so to investigate the origin of the low gas masses in red 
early- type galaxies. In light of this, we now consider the 
results of our observations in both the stacked and indi- 
vidual analysis. We show the results of our observations 
in Figure m In each case, we plot the mass of molecular 
gas as a function of the time elapsed since the start of 
the last significant episode of star formation ty from S07. 

In the case of non-detections, we must make an as- 
sumption on the expected line profile to determine the 
amount of gas that would lead to a detection. We divide 
our analysis into two separate avenues here: first, we per- 
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Fig. 2. — CO(l — » 0) line profiles for our targets. We plot the temperature as a function of velocity measured from the redshift of the host 
galaxy. For each object, we give the SDSS name and emission line classification (SF, AGN+SF, Seyfert). In the top four rows, we show the 
detections, while in the bottom row, we show a sample of three non-detections. For the detections, we indicate the best-fit line profile from 
our CLASS data reduction. Note that all detected lines are relatively narrow (~ 70 kms~^) and many objects show double-peaked profiles. 



form a very conservative analysis using stacked spectra to 
establish that we see a trend between age/emission line 
class and gas fraction. After establishing that a trend ex- 
ists, we then consider individual galaxies using less con- 
servative upper limits. 

4.1. Analysis of stacked spectra 

We wish to establish the most conservative possible 
upper limits on our non-detections. Despite the low rms 
noise levels of our spectra, it is conceivable that a very 
broad line might lurk just below our sensitivity limits. 
In order to guard ourselves against this, we stack the 
non-detected spectra in each emission line class to in- 
crease our signal-to-noise and assume the broadest line 



that can reasonably be expected, a line with FWHM of 
300 kms""'^. From this, we obtain 2a gas mass upper 
limits for the non-detected sources in each class of 10.0, 
3.4 and 1.7 x 10^ Mq for the SF, AGN+SF and Seyfert 
classes. The stacked spectra for the SF, AGN-I-SF and 
Seyferts are shown in Figure [31 they show no evidence 
for any lines. For the objects with detected lines, we sim- 
ply take the mean of the derived gas masses, as stacking 
might distort the line profile and thus the derived masses. 

In panel (a) of Figure IH we show the mean gas masses 
of those galaxies which are detected in both the SF and 
the AGN-I-SF classes. We furthermore show the upper 
limits for the stacked non-detections assuming a very 
conservative line width of 300 kms^^. The starforming 



Molecular Gas Destruction by AGN Feedback 



5 



TABLE 1 
SDSS & IRAM Data 



SDSS Object Name R.A. Doc. Emission Stellar Mass Total CO Flux^ Molecular Gas Age of 

Line Classl M,t„i,„ ScO(l-.0) ""^^^ '''s^" Starburst* ty 

(J2000) (.12000) (XIO^" Mq) (K kms^l) (xlo'* Mq ) (Myr) 

^ TT, 



SDSS 


J075 109 


77+342636 




















SF 






+ 


















180™ 
^8U_7Q 




J u -1- .^uoy 


73 1 142S24 












+ 14 














+ 






63+019 






+ 






70*0 


SDSS 


J002227 














+ 13 








SF 






+ 


















340™ 




J U -1- -1- 4*00 


^ _L 1 AOTOfi 
00-f--l-'**i 1 












+ 14 








AG 1*1 1 £F 






+ 






















J UUUOOD 














+ 14 








AGN + SF 






+ 


















iqqnl'lSO 
199U_^720 




J 100124 














+ 37 








AG 1*1 1 £F 






+ 






28 i 055 






+ 






460200 
4DO_2oo 


SDSS 




ux-poD*iy*iij 












+ 36 














+ 






73 + 10 






+ 






200™7o 




Jl 10234 


riT 1 .1 Q n n 
y * -p^jSjiOUy 












+ 42 








AGn 1 £F 






+ 






56 + 11 






+ 






1 50*0 




mani o'^ 
J uy u -1- ^0 




























+ 






39 i 070 






+ 






1 90*0 




J Uo^ r jsy 


-t n i -jci-ii-ioo 
^y-p ODUU^^ 












+ 36 








AGN4-SF 






+ 






97 + 16 






+ 






1 50IOO 
1--^0_9Q 


SDSS 


J082104 


65 + 412108 


5 


08 


21 


04 


7 


+ 41 


21 


OS 


6 


AGN + SF 


3 


5 


± 





63 


< 0.07 








1 


65 


3I0IOO „ 
^^—130 


SDSS 


,7102847 


41 + 404231 


3 


10 


28 


47 


4 


+ 40 


42 


31 


4 


SF 


3 


8 


± 


1 


11 


0.53 + 0.1 




3 


+ 


1 


00 


1 90*0 
^^0_4Q 


SDSS 


J124450 


01 + 422123 


3 


12 


44 


50 





+ 42 


21 


23 


4 


SF 


8 


3 


± 


1 


69 


1.38 ± 0.190 


13 


4 


+ 


1 


84 


I40IOO 

^^0_gQ 


SDSS 


J125716 


71 + 424625 


9 


12 


57 


16 


7 


+ 42 


46 


26 





SF 


3 


1 


± 





19 


< 0.21 






< 


4 


09 


,,nl20 


SDSS 


J075025 


78+444621 


3 


07 


50 


25 


8 


+ 44 


46 


21 


3 


Seyfcrt 


2 


8 


± 





31 


< 0.07 






< 


1 


18 


47nis>0 

*^0_i4o 


SDSS 


J111850 


04+422541 


8 


11 


18 


50 
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25 


41 


8 


SF 


2 


6 
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47 


0.55 + 0.10 


6 


2 


+ 


1 


13 


17o5'0 


SDSS 
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81 + 364803 


7 


08 


40 


14 


8 
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48 


03 


7 
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3 


1 
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58 
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< 


1 


45 
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SDSS 
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5 


10 


38 


55 


9 
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21 


57 


6 


Seyfcrt 


1 


7 


± 





24 


< 0.07 






< 


1 


40 


OOoi^gOgo 


SDSS 


J093014 


24+360909 


5 


09 


30 


14 


2 


+ 36 


09 


09 


6 


AGN + SF 


2 


6 


± 





53 


0.17 ± 0.065 


2 


1 


+ 





79 


100™gg 


SDSS 


J082115 


75 + 355924 


6 


08 


21 


15 


8 


+ 35 


59 


24 


7 




4 





± 





41 


< 0.07 






< 


1 


68 




SDSS 


J082507 


10+410319 


4 


08 


25 


07 


1 


+ 41 


03 


19 


5 


AGN + SF 


3 


8 


± 





60 


< 0.07 






< 


1 


51 


220l''0 
■'''"-120 


SDSS 


J084216 


23+360141 


4 


08 


42 


16 


2 


+ 36 


01 


41 


5 


Seyfcrt 


3 


8 


± 





35 


< 0.07 






< 


1 


37 


1740l_020^„ 


SDSS 


J104705 


19+411514 


5 


10 


47 


05 


2 


+ 41 


15 


14 


5 


SF 


1 


7 


± 





31 


0.34 ± 0.090 


3 


4 


+ 





89 


190«40 


SDSS 


,7103843 


56 + 404221 





10 


38 


43 


6 


+ 40 


42 


21 





AGN + SF 


2 


8 


± 





49 


< 0.07 






< 


1 


39 


8903^0 



1 BPT emission line classification, described in detail in Schawinski et al. (2007). 

^ Upper limits for the total CO flux are based on the rms noise of the smoothed spectrum and assuming a line of width 70 kms^^, they arc Icr limits. 

^ Molecular gas masses are derived from the total CO luminosity and assuming a conversion from CO to molecular gas mass of a = 1.5 K km s^^ pc^ (Evans ct al. 

2005). The gas mass upper limits are 2ct as in the Figures. 

^ The age of the last episode of star formation ty is taken from Schawinski ct al. (2007) and is based on fitting the U V/optical/NIR spectral energy distribution and the 
stellar (Lick) absorption indices. 



TABLE 2 
IRAM Data Line Properties 



SDSS Object Name-^ 




Line Width 
kms-1 




Li 




Peak 
K 


SDSS 


,1012039 


73+142824 


2 


120 + 20 


5 


1 




10 


-i 


SDSS 


.1100124 


35 + 375046 


8 


60 ± 25 


2 


4 


X 


10 


-3 










35 ± 12 


3 


3 


X 


10 


-3 


SDSS 


.7083539 


01 + 362923 


2 


55 ± 18 


4 


5 


X 


10 
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5 


1 


X 


10 


-3 


SDSS 
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2 
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3 


2 


X 


10 
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71 ± 28 


2 


8 


X 


10 


-3 


SDSS 
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4 


5 


X 


10 


-3 
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1 
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early-type galaxies at young ages have high molecular gas 
masses of the order of up to ^ 10^ Mq. On the other 
hand, the pure (Seyfert) AGN where no traces of ongoing 
star formation are left in the optical emission hues are 
undetected in CO with upper limits on the gas mass sig- 
nificantly lower than the younger SF objects. In between 
lie the AGN-I-SF objects: the detections are systemati- 
cally younger and have mean gas masses slightly below 
the SF objects. The non-detected AGN+SF objects are 
older than their detected counterparts, but younger than 
the Seyferts. During in the AGN+SF phase, at about 200 



Myr after the onset of star formation, the molecular gas 
mass drops precipitously. 

It should be noted that the presence of CO correlates 
with other measures of star formation and so galax- 
ies with more CO have more vigorous star formation 
(lKennic utt"l998t lGao fc Solomonll2004l;lNaravanan et al.l 
120051) . However, the star formation history parameters 
derived by S07 would account for the implied larger mass 
fraction in the young population and still return the cor- 
rect age for it. 

4.2. Analysis of individual spectra 

On the basis that we can establish the significance 
of the observed trend in panel (a) of Figure [5] using 
the stacked spectra, we now turn to invididual galax- 
ies. Since S07 resolve the age of the current episode of 
star formation and provide emission line classes for each 
individual object, we can also plot these to reveal trends 
within populations. Half of our detections shows sin- 
gle, narrow C0(1 0) lines. The other half shows two 
components. The median FWHM of the single lines and 
the individual components is 70 kms~^. These relatively 
moderate widths are consistent with the fact that the 
galaxies in our sample are only intermediate-mass ob- 
jects and that they are undergoing minor residual star 
formation. They are not undergoing major starbursts 
requiring enormous gas reservoirs, nor are they spiral 
galaxies with extended disks with broad double-horned 
lines. 

In the absence of comparable samples, we use the me- 
dian FWHM of our sources with detections of 70 kms~^ 
to estimate the upper limts on the gas mass in the non- 
detections. These are not as rigorous as the limits of the 
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Fig. 3. — CO(l — > 0) stacked spectra for the Seyferts and the 
AGN+SF non-detections. Assuming a very broad line as a worst- 
case scenario with a FWHM of 300 kms~^, the 2cr gas mass upper 
limit of the SF, SF-f AGN and Seyfert stacks are 10.0, 3.4 and 
1.7 X 10* Mq. 



stacked spectra, but as we will show, the results are con- 
sistent. With this analysis, we neglect the possibility of 
double-horned profiles in the analysis, which would in- 
crease the upper limits of the non-detections by up to 
a factor 2. The analysis of the stacked spectra however 
does not support this latter scenario. 

The individual galaxies show the same trend and al- 
low us some insight into variations within each popula- 
tion. We present these in panel (b) of Figure [4] and the 
gas masses and ages ty are given in Table [TJ We see 
here that the SF molecular gas masses are within a fac- 
tor of a few of each other and reach up to ^ 10^ Mq. 
All three Seyfert AGN are individually undetected. In 
between are the AGN-I-SF objects, two of which have 
high masses with the remainder either undetected or sig- 
nificantly lower than the starforming galaxies. Figure 
shows the same trend as Figure 2^. The individual 
gas masses and upper limits further support the observa- 
tion that the gas mass drops precipitously at ^ 200 Myr 



within the AGN-f SF phase. 

To assess the significance of the difference between the 
populations we divide the sample into two groups, one 
younger than 200 Myr and one older. We use Bayes' 
theorem to find the posterior distribution of the pop- 
ulation mean and variance of the two groups' molecu- 
lar gas masses by assuming that each group is drawn 
from a different normal distribution. Samples on which 
only upper limits were found are described as gaussian 
data points with zero mean and the appropriate width. 
We place a fiat prior on the population means and stan- 
dard deviations. After computing the joint posteriors, 
we marginalize over the standard deviation. The result 
is a probability distribution for each of the two popula- 
tions means regardless of their variances. Together the 
distributions show that the younger group has a higher 
mean with > 95% probability. 

4.3. The case for AGN feedback: the Schmidt Law and 
gas depletion time scales 

Such a sudden drop in molecular gas mass as we ob- 
serve cannot be accounted for by star formation alone. 
To demonstrate this, we compute the evolution of the 
molecular gas mass for both the mean SF galaxy and 
each of the SF early-types with C0(1 0) detections , 
using the Schmi dt Law for star formation (jSchmidtill959l : 
lKeniiicuttlll998f) : 



SFR 



gas 



(1) 



We assume a fiducial efficiency of 2% (e = 0.02), which 
includes the effect of SN feedback. This law appears 
to be universally applicable a cross galaxy mo rphologies, 
including early- type galaxies (jCombes et al.|[2007D . We 
compute the dynamical times of each galaxy from its 
radius and stellar mass (tdyn — \/ B? I2GM)\ for the 
stacked SF galaxies, we use the median tdyn of 50 Myr 
of the SF galaxies. We show the resulting evolutionary 
tracks in Figure [4^ and c (dashed lines). 

The depletion time scales derived this way from the 
Schmidt Law are on the order of several Gyrs, much 
longer than the observed timescale for gas depletion. The 
tracks shown in Figure [4^ and c are inconsistent with 
a scenario where star formation alone depletes the gas 
reservoir. 

Would it be possible to adjust the parameters used to 
result in tracks for star formation alone that are consis- 
tent with our observations? Adjusting the mass and ra- 
dius to achieve the rapid gas reservoir exhaustion we ob- 
serve results in parameters not characteristic of a galaxy- 
wide starburst, but perhaps a circumnuclear ring. Al- 
ternatively, a drastic increase of the star for mation ef- 
ficiency e towards unity driven by an AGN (lsiia[2005h 
might suffice, though this completes the circle and brings 
us back to AGN feedback. Our modeling does not ac- 
count for the amount of recycled gas from stellar mass 
loss or further gas accreted, both of which would further 
increase the amount of molecular gas available over time 
and thus extend the depletion time scales significantly. 
This caveat further strengthes our argument against star 
formation alone. 

This decrease in molecular gas coincides with the time 
when galaxies transition from pure star formation to an 
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Fig. 4. — The observ ed molecular gas masses Mhj as a function of ty, the age of the last significant episode of star formation 
llSchawinski et ani2007bh . The blue stars are early- type galaxies classified as starforming, the green squares are AGN+SF composites 
and the red circles are Seyfert AGN. In the top panel, we show the results from the analysis of the stacked spectra. The upper limits are 
2rT assuming a typical linewidth of 300 kms~^. In the bottom panels, (b) and (c), we show the individual gas masses for each galaxy. The 
upper limits are 2a assuming a typical linewidth of 70 kms~^. The values of ty are based on assuming an e-folding time of 100 Myr for 
the starburst; allowing this time to vary simply extends the error bars, but does not move the best-fit values. Instantaneous bursts or 
constant star formation models yield poor fits to the data. In panels (a) and (c), we show the tracks showing the future evolution of the 
gas mass assuming the Schmidt Law for star formation. In (a), we show this track for the mean SF galaxy, wherease in (c), we show tracks 
for each individually detected SF galaxy. These tracks illustrate that given an evolution driven purely by star formation, the depletion of 
the gas reservoirs on the short time scales we observe is not feasible. These tracks assumes a dynamical time scale of r = 100 Myr and 
is inconsistent with the systematically lower gas masses of the AGN-fSF objects — whether detected or not. We mark the point in time 
ty = 200 Myr when the low-luminosity AGN in the AGN+SF objects typical switch on. 



AGN-I-SF composite as the AGN switches on of ~ 200 
Myr (SOT). Figure [5] shows this same trend in the emis- 
sion line diagnostic diagrams. There the two AGN-I-SF 
composites with the gas masses similar to the starforming 
objects reside in close proximity to them; those AGN-f-SF 
objects further away from the starforming locus are ei- 
ther undetected or have significantly lower gas masses. 

5. DISCUSSION 

Models of hierarchical galaxy formation suggest that 
AGN feedback plays an important role in the forma- 
tion and late-time evolution of early-type galaxies. Cur- 
rent simulations invoke two modes of bla ck hole growth 
and AGN feedback: a 'quasar mode' (e.g. iSpringel et al.l 
I2005t iDi Matteo et all 120051: iHop kins et al."2008') and a 
'radio mode' (e.g. ICroton et al. [ [2006; .Bower et al...2006. : 
ISiiacki et al.][2007l) . 



5.1. Two modes of AGN feedback 

The qua sar mode, imp l ement ed in simulations such 
as that of ISpringel et al.l (|2005f ). assumes that during 
a major merger a fraction of the rest mass of the gas 
accreted onto the central black hole is liberated as en- 
ergy injected back isotropically into the gas reservoir 
of the host galax y. Simulations of gas-rich mergers at 
very high redshift (iNaravanan et al.l l2006l: iLi et al.|[2007t 
[Narayanan et al.l l2008f l indicate that quasars might be 
driving massive, powerful outflows as the energy lib- 
erated by black hole accretion is coupled to the gas 
in the form of thermal energy. These outflows may 
succeed in halting s t ar fo rmation on short time scales 
(|Hamann fc Ferlandl |1992| ). as required for the most 
massive early- type galaxies ([Thomas et al.l I2005D . On 
the other hand, observations of quasars at high red- 
shift yield very high CO luminosities and hence high in- 




Fig. 5. — The objects in our sample on emission line diagrams lIBaldwin et al.l 119811 : IKewlev et al.l 120011: IK auffmann et al.l 120031 : 
IKewlev et al. 1'2006'). In panel (a), we delineate starforming galaxies by the dashed line fKa03: iKau^nann et aTn2Q03). B etween the 
dashed and solid line (KeOl) reside AGN+SF composites, while beyond the solid line there are the pure AGN (IKewlev et al.l l266iV The 
straight line (S07) divides Seyfert AGN from LINERs. The blue stars represent starforming early-type galaxies, the green squares AGN+SF 
composites and the red circles Seyfert AGN. Filled symbols are CO(l — > 0) detections and open symbols are non-detections. The size of 
the symbol scales with the logarithm of the molecular gas mass (see legend in (a)). In panels (b) and (c), we show a similar line ratio 
diagram using the classification of Kcwlcy ct al. (2006). The species [SII] and [OI] used in (b) and (c) are sensitive to low ionisation states 
and complement (a) based on [Nil]. Note that the two AGN-I-SF galaxies with large molecular gas masses are also modest radio sources. 



ferred gas masses (|Cox et al.ll2002l: iBertoldi et al.l 120031: 
I Walter et al.ll2004 ). which suggests that feedback in these 
objects does not necessarily result in the instantaneous 
removal of molecular gas fr om the s ystem . 

The r adio mode (Crot on et al.1 ()2006l ): iBower et al.1 
(|2006( ): iKhalatvan et all (|2'00l ). instead, is important 
for the evolution of galaxies and is supposed to occur at 
later times. A low-level AGN is thought to be injecting 
a fraction of the rest mass energy of the hot gas accreted 
by the central black hole into the ISM over an extended 
period of time to keep the host galaxy quiescent. 

5.2. Molecular gas reservoir destruction by 
low-luminosity AGN 

Our galaxies help us illuminate the role of AGN in 
suppressing late-time star formation. In previous work 
(S07) we have identified a sample of low- and interme- 
diate mass early-type galaxies that display an evolution- 
ary sequence from minor star formation through AGN to 
the red sequence, most likely driven by AGN feedback. 
In the present paper we investigate through IRAM mm- 
observations the presence and absence of cold, molecular 
gas reservoirs in objects selected along this sequence. 

The results presented here support the presence of late- 
time AGN feedback invoked in the simulations. Our ob- 
servations suggest a scenario in which molecular gas is 
removed or heated ~ 200 Myr after the peak of recent 
star formation in our systems. Most importantly, this ob- 
served drop coincides with the timescale on which these 
galaxies transition from being primarily star forming to 
being dominated by AGN. We show that significantly 
longer timescales would be expected if the molecular gas 
was to be consumed by star formation alone. 

The AGN detected in our sample are low-luminosity 
AGN. They are neither powerful quasars nor radio- 
loud. Their accretion rates and efhciencies are low 
(106 < L[o„i] < 108 Lq, -3 < log(L[oiii]/fT^) < 



0; S07)^°. Out of our entire sample only two tar- 
gets are detecte d at 1.4 GHz by t he Very Large Array 
FIRST Survey (jBecker et al.lll995[ ) at a detection limit 
of 3 X 10^1 W Hz-i, both of which are AGN-hSF ob- 
jects with modest L1.4 ghz of 6.7 x 10^^ W Hz"^ and 
2.5 X 10^2 W Hz-^ These two AGN-I-SF objects are 
those with the two highest molecular gas masses. 

Some theoretical models have discussed scenarios 
where low-luminosity, low accretion efficiency AGN 
of the kind we see here regulate star formatio n 
()Hopkins fc HernauistI 120061: ICiotti fc Ostrikeil [2007'). 
Furthermore, the time delay between peak SF and AGN 
Eddington ratio and the moment of 'blowout' - when the 
AGN drives out the gas reser v oir - i s present i n the mod- 
els of iHopkins &: Hernguisti (|2006D and Ho pkins et al.l 
(|2006f) . The process we see thus might be physically sim- 
ilar to quasar mode feedback as envisioned by theorists, 
but fulfills the task of the radio mode in moderate-mass 
early-type galaxies of suppressing further late-time star 
formation. 

We conclude that low-level AGN appear to be powerful 
enough to destroy the molecular gas reservoir and sup- 
press star formation in early-type galaxies. Furthermore, 
we note that along the time sequence established and de- 
scribed by S07, the peak of the AGN accretion rate and 
efficiency occur during the Seyfert phase; the AGN phase 
that suppresses star formation is less efficient than this 
peak. 

5.3. Suppression vs Truncation 

At high redshift, early- type galaxies (or their progen- 
itors) experience vigorous star formation. The large en- 
ergy output from quasar-like AGN is required to bring 

The quantity iioni]/""*: introduced bv IKewlev et all (120061 ). 
is a proxy of the Eddington ratio. L[oiu] traces the accretion 
rate, while cr^ traces the black hole mass via the M, — cr relation 
(IGebhardt et al.ll2000l : IFerrarese fc MerritUl2000h . 
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such massive star formation to a halt. We refer to this 
type of AGN feedback as 'truncation mode'. In order to 
account for the observed high star formation efficiency 
in these galaxies, it is possible that a 'trigger mode' may 
precede the truncation mode (Silk 2005). 

Clearly, the process we are probing with the present 
work is very different. We see AGN feedback operat- 
ing at late times in early-type galaxies associated with 
low-luminosity AGN. No quasar-like activity or powerful 
radio jets are observed. We witness the suppression of 
star formation in early- type galaxies at recent epochs. 
We therefore prefer to call this type of AGN feedback 
'suppression mode' as opposed to the truncation mode 
at early times. The low- and intermediate mass galax- 
ies studied in our undergo very minor episodes of star 
formation and thus allow us to study this suppression 
process. Apparently, the AGN in this lower mass range 
is not powerful eno ugh to inhibit residual s tar formation 
completely (see also lSchawinski et al.ll2006l ). As shown in 
Figure m however, without the energy input from AGN, 
those episodes of star formation would be more signifi- 
cant and would drive the objects even further away from 
the red sequence. The observed suppression by AGN 
feedback ensures that these objects will join or return 
to passive evolution on the red sequence. Note that the 
massive early-type galaxies in our sample, instead, do 
not show any signs of residual star formation. As dis- 
cussed in S07, the LINER activity detected in some of 
them might be the 'smoking gun' of highly effective sup- 
pression mode AGN feedback. 

Finally, it is interesting to note that the suppression 
mode detected here acts on timescales around several 
hundred Myr, which is the typical timescale of an AGN 
duty cycle. Hence, suppression of residual star formation 
in early-type galaxies at late times ma y not be a contin- 
uous but a periodic process (see e.g. ICiotti fc Ostrikeri 
l2nf)l . 

6. CONCLUSIONS 

The molecular gas reservoirs fuelling late time star for- 
mation in early-type galaxies must be destroyed to ac- 
count for their passive evolution on the red sequence. 
This suppression of star formation in early-type galaxies 
is now commonly attributed to AGN feedback wherein 
the reservoir of gas is heated and expelled by the en- 
ergy input from an accretion phase onto the central su- 
permassive black hole. In this Paper, we present ob- 
servational evidence for this process occurring in low- 
to intermediate-mass early-type galaxies at low redshift 
and identify low-luminosity AGN - not quasars or radio 
galaxies - as the culprit. The systems studied here are 
not the massive galaxies in clusters where the radio mode 
is generally invoked to suppress cooling. 

In order to empirically determine whether and at which 
stage AGN are involved in the destruction of significant 
molecular gas reservoirs, we perform an observational in- 
vestigation. We have observed a sample of low-redshift 
intermediate mass early-type galaxies during the process 
of moving from the blue cloud to the red sequence via a 
low-luminosity AGN phase. We observed galaxies along 
the time sequence of S07 with the IRAM 30m telescope 
near Granada, Spain to determine the amounts of molec- 
ular gas present via the C0(1 — > 0) line. 

We perform two separate analyses: First, we take 



the average gas masses in each of the SF, SF-fAGN 
and Seyferts galaxies and stack the spectra of the non- 
detections to increase signal-to-noise and therefore the 
upper limit on the gas mass. For the non-detections, 
we assume an extremely conservative line width of 
300 kms~^ and find a statistically significant trend. The 
gas mass in the SF galaxies drops approximately 200 Myr 
after the start of star formation, during the AGN-I-SF 
phase, and is far too rapid to be accounted for by star 
formation alone. None of the four Seyfert AGN were de- 
tected, suggesting that the gas reservoirs were already de- 
stroyed in the preceding AGN-I-SF phase, i.e. prior to the 
peak AGN luminosity and accretion efficiency (c.f. S07). 
The molecular gas reservoir drops by at least an order of 
magnitude in mass in a very short time. The Schmidt 
law for star formation implies depletion timescales of sev- 
eral Gyr, even if contributions such as further accretion, 
cooling and mass-loss by evolved stars is ignored. We 
thus interpret the low-luminosity AGN in the AGN+SF 
phase as the culprit responsible for the destruction of 
the gas reservoir. It is possible that the feedback work of 
the AGN might have begun in an obscured fashion dur- 
ing the SF phase, during which the feedback from stars 
might have aided the gas removal. 

Second, we perform the same analysis for each galaxy 
assuming a measured line width of 70 kms~^ for the non- 
detected galaxies. We find the same trend: starforming 
early-type galaxies have substantial molecular gas reser- 
voirs. A few AGN-I-SF early-types are detected with 
diminished gas reservoirs, while others are undetected, 
with gas masses inferred to be at least an order of mag- 
nitude lower. The individual galaxy data points show 
that the depletion of the molecular gas reservoir occurs 
within the AGN-I-SF phase about 200 Myr after the start 
of star formation. The Seyfert AGN are all undetected. 
For the individual galaxies, we perform a Bayesian signif- 
icance test and find significant evidence for populations 
of galaxies with star formation more recent than 200 Myr 
and older than this value having different molecular gas 
masses. 

We then discuss the implications of our result and how 
it impacts our understanding of AGN feedback. The pro- 
cess we see occurring suppresses residual star formation 
in intermediate-mass, low redshift galaxies and so could 
be termed a 'suppression mode' of the AGN, which main- 
tains the red colors of early-type galaxies by suppressing 
residual star formation. This mode stands in contrast 
to the 'truncation mode' which is thought to shut down 
star formation during the formation of massive galaxies 
at high redshift. 

The AGN phase responsible for this suppression mode 
is a low-luminosity event and not a powerful quasar 
or radio galaxy phase. The fact that we observe low- 
luminosity AGN impacting the fuel for star formation in 
early-type galaxies is an important step towards under- 
standing galaxy formation. 
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